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SUMMARY 



Problems of hydraulic flow resistance and heat 
transfer for streams with velocities comparable with 
acoustic have present great importance for various ■' 
fields of technical science. Especially, they have 
great importance for the field of heat transfer in de- 
signing and constructing boilers-of the "Velox" type. 
In this article a description of experiments and their 
results as 'regards definition of the laws of heat trans- 
fer in differential form for high velocity air streams 
inside smooth tubes are given, 

EXPERIMENTAL APPARATUS AND PROCEDURE 
-' General "Scheme of Experimental Apparatus 
and Experimental Section 



The general scheme of experimental apparatus is 
shown in figure 1, The air came in from the compressor 
into the apparatus through the large surge tank B, 
which has the purpose of. removing the pulsations and 
separating the water fog. ' The orifice D was installed 
to' measure the air flow rate and the air passing through 
was directed into' a zig zag conduit ,P inside of which 
there was an electric heater. The heated air passed 
through the mixing section of the pipe G into the ex- 
perimental tube A x where it lost a part of its heat to 

ac -. - — — .-■ — .. T . - , . r - . , - ■ ■ i . . | - |- | n | _ 

*Jour. Tech, Phys. (USSR), vol. -IX, no. 9, .1939, pp. 808-818. 

NOTE: Translation received from Massachusetts Institute of 
Technology, Cambridge, Mass*- 
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cooling water which flowed' through' the casing around the 
experimental tube, losing its pressure and expanding up 
to the limit of atmospheric exhaust pressure. The air 
discharged from the experimental section at high velocity 
which approached acoustic velocity. 

The experimental section is shown in figure 2. A 
smooth copper tube of 14— millimeter inside diameter and 
16— millime ter outside diameter "by 1431— m illimeter length 
was placed axially inside a brass tube 19 . 8— mill ime t er 
inside diameter. The heated air flowed with high veloc- 
ity through the inside tube and the cooling water flowed 
in the annular space. Transition from the large pipe to 
the exper iir.ental tube was accomplished smoothly, by means 
of a slowly converging nozzle, A length of 21 millimeters 
at the exhaust end of the experimental tube was given a 
diverging oonic shape. Thus the experimental tube may be 
considered as a nozzle, first converging and then after 
a long narrow throat, enlarging again. Tlow velocities 
as high as acoustic velocity and corresponding pressure 
relations were achieved at the point of transition from 
the cylindrical section to the conic section. By this 
means the opportunity to make relatively easy all measure- 
ments, including even those in the acoustic region were 
obtained. The air pressure was measured at a number of 
points along the experimental section. The temperatures 
of the inside wall of the air tube as well as those of the 
cooling wat-e : r, we r,e also measured along the flow. In the 
experimental tube 9 static pressure tops of 0. 08— millimeter 
diameter were located to measure the changes in static 
pressure as the fluid flowed through the tube. In order 
to eliminate roughness formed by boring the pressure taps 
a special drill by which the -burrs were cut off from in- 
side the tube was made. 

After- the holes were. made and the pressure taps - 
attacked the transverse section of the tube and its vari- 
ation all along 1 the length-was measured by gradually fill- 
ing- the tube with mercury in the section between two pres- 
sure taps and weighing the -mercury. -The measurements' 
showed that the diameter of the tube was sufficiently 
uniform all along. The average , inside area of a trans- 
verse section of the tube was 155,3 square millimeters 
and the mean diameter 14.07 millimeters. At seven points 
in - the outside wall' of the copper tube thermocouples were 
inserted to measure the tube— wal-1 temperatures and. their 
variation along the tube. t J!or this purpose, small grooves 
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„,w.e t re made in the tube Wall,' the'-'two thermocouple wires 
were inserted in the gr o oVe'V and ' soi'd'e red • over flush 
with the tube wall.' "In order to detect any asymmetry 
of cooling of the tube which might have occurred during 
the experiment, four thermocouples were inserted along a 
line parallel to the tube axis and th'ree alternate couples 
were placed along a line diametrically opposite. The dis- 
tances Between the thermocouples can. he seen in figure 1. 
To measure the temperature of the cooling water at points 
along the flow inside, thermocouples of wire of 0. 2-milli- 
meter diameters were placed in the annular space. The 
diameter of the smoothed soldered junctions was 1.0 milli- 
meter.. Each thermocouple was tied by a thin wire close 
to the' experimental tube so that the free end of the couple 
was slightly above the tube surface. A thermocouple in- 
stalled in this manner with relatively large thickness of 
junction equal to one— half the crack width assured measure- 
ment of the bulk temperature of the water flowing at this 
point. There were six such thermocouples along the tube 
in addition to two others to measure the inlet and exit 
temperatures of the water. These thermocouples were dis- 
tributed similarly to the tube— wall thermocouples. All 
the thermocouples for measurement of the wall and water 
temperatures were carefully insulated by a special varnish 
and calibrated. The water thermocouple distribution in 
the annular space is shown in figure 1, 

In making up the experimental section it was necess-a— 
ry to keep the experimental tube centered inside the 
casing to insure symmetrical flow of cooling medium. This 
was done by using, as spacers, three small wires 20 -milli- 
meters long and of diameter equal to the crack width sol- 
dered to the tube' wall at each of these points. Besides, 
this, centering rings were soldered on each end of the ex- 
perimental tube and the ends of the casing were then sol- 
dered to these rings. In order to avoid any deformation of 
the experimental section because of the differences in ther— a I 
mal expansion of experimental tube and casing, the latter 
was cut transversely in the t middle into two sections. The 
ends of the casing could then move freely inside a metal 
ring placed over the cut in the casing. Tightening of the 
joint was accomplished, by use of a rubber tube placed over 
the ring and wired to both ends of the cut casing;' thus 
the rubber served to absorb all the deformation resulting 
from thermal expansio*n of tube and casing. The thermo- 
couple leads were brought out through small holes in the 
casing. Short tubes of 3 millimeters out s ide— diame t er 
were carefully soldered to the pressure taps in the test 
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section. These tubes were brpught .out through holes made 
in the casing as shown in figure 2 "by -shorty thick— walled, 
metal— armored rubber tubes placed over the small tubes. 
A tube passing through and soldered to the casing was in- 
serted into the outer end of the rubber. Thus there was 
assured the necessary compactness and elasticity of the 
joints to the manometer and all 'possible security in 
measurement of the static pressure of the air. To the 
outside wall of the casing there were also soldered five 
thermocouples to measure the casing temperature along the 
experimental section. The casing was insulated with as- 
bestos to avoid heat loss. 



.Experimental Procedure and Measurements 

Before the experiment was begun the air and the cool- 
ing water were started flowing in the apparatus and the 
electric heater was turned on to heat the whole apparatus 
and to reach steady—state operation. After the desired 
flow conditions were reached, the experiment was started. 
Before and during the experiment steady. air flow rates 
were maintained by continuous observation of the level 
of the differential water manometer on the orifice and by 
careful regulation of the valve B. The duration of one 
experiment varied from 15 to 30 minutes. During this 
time all data were taken three times. Some of the exper- 
iments were repeated from two to three times under the 
same conditions. Because the air flow rates were regu- 
lated, the pressure in the experimental tube varied in- 
significantly during the experiment. 

In order to assure even distribution of cooling water 
around the annular space, aritificial turbulence was intro 
duced into the water flow by admixture of compressed air 
which was put into the water line before the water entered 
the casing. The well— mixed air— water emulsion passed 
through the annular space at high velocity and thus during 
the experiments uniform smooth curves of thermocouple read 
ings were obtained for the cooling water and the experimen 
tal tube. The inevitable inaccuracy in heat quantities 
due to the air in the water is insignificant even if the 
volume of the air were three times that actually used. In 
our experiments c ount er— curr en t flow of air and cooling 
emulsionwere used. With parallel flow the main portion of 
the change in water flow would occur in the first portion 
of the experimental section and this first section inter- 
ests us least because it represents only the beginning of 
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the velocity range o-f-. int-e rest , With counter— flow the 
change in water temperature is more' uniform with insig- 
nificant con-cavity toward the axis of length. , In this 
latter case the curve of water temperature from inlet 
to exit of the casing is almost linear along the entire 
experimental section. This assured greater accuracy in 
differentiation of these curves. 

The distribution of air pressure along the flow was 
measured "by a mercury differential manometer having 10 
legs, A detailed description of this manometer is given 
in the work on measurement of hydraulic resistance. Nine 
legs of the manometer were attached to the pressure taps 
and the tenth point communicated to the atmosphere. Thus 
the differential manometer gave a distribution of mercury 
levels corresponding to the pressure gradient in the tube 
above atmospheric pressure. The air flow was measured 
by a standard orifice made of stainless steel placed in 
a 70— millimeter ins ide— diameter pipe upstream from the 
electric heater. The orifice diameter was 21 millimeters. 
The pressure drop across the orifice was measured by a 
U— tube manometer filled with colored water. The absolute 
pressure upstream from the orifice was measured by a pre- 
cision Bourdon gage to six atmospheres with a scale grad- 
uated to 0.03 atmosphere. The air temperature upstream 
from the orifice was measured by a thermocouple. 

The quantity of heat Q, (per unit time) lost by the 
air in the test section was determined by measuring the 
weight rate of flow of cooling water and temperature rise 
of the water. In addition, the curve of water temperature 
against length of test section was measured so that by 
differentiation of this curve the instantaneous rate of 
heat transfer could be determined. 



dL 



cG- at 
dL 



(1) 



where 



G e mass flow rate of water, kg/hr 

dt 6 /dL* derivative of water temperature t 6 with respect 
to length, deg.C/meter 



*1 was changed to L in translation to facilitate reading 
the typewritten text. 
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dCj/dL 



derivative of with respect to length in 

cal/(m)(hr) 

specific heat of water = 1 cal/(g)(°C) 



c 



The entire heat loss from the air along the test 
section was also computed from measurement of air flow 
rates and inlet and exit temperature and pressure for 
the test section. With these mass flow rales and the 
same inlet and exit pressures and temperatures obtained 
in the 3— inch sections preceding and following the test 
section, the heat loss from the air was measured in a 
3— inch ins ide— d iamet er pipe (where changes in kinetic 
energy of the air were negligible ) . The two independent 
measurements of Q, checked in our experiments to with- 
in 5 percent. 



In all our experiments at the exhaust end of the 
tube acoustic velocity at that temperature >as attained, 
.Table I shows the complete measured data fcr all experi- 
ments. Prom these data plots of measured quantities were 
made. 

Figure 3 shows curves of the pressure distribution. 
The abscissa is marked as distance in millimeters from the 
entrance of the test section and the ordinate as the pres- 
sure in millimeters of water, .trur each experiment plots 
of temperatures of inside wnlL of test section, cooling 
water temperaturej and casing wal?_ temperature have been 
made and these are shown in figure 4, 

The motion of a compressible gas is described by 
the following formulas: 

( 1) Conservation of matter 



EXPERIMENTAL DATA AND THEIR MEANING 



G = 



Pw 



(2) 



V 



(2) Equation of state 



pV » RT 



(3) 
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(z) Energy "balance 

1*1 + 6p I + Q = ?p In + A ^2- (4) 

where 

G air flow rate, mass per unit time 

I* transverse area of tube 

w flow velocity 

V specific volume 

p absolute pressure 

I absolute temperature 

w, V, P, T mean values over the section* 

Cp ' specific heat at constant pressure 

B, gas constant 

A thermal equivalent of mechanical work 

g acceleration of gravity 

Q, heat output per kilogram of air in the test 

section from the inlet to the section in 
quest ion 

Sub script: 

o inlet section 

If started at the large tube "before the test section 
the initial kinetic energy A <-w Q 2 /2g) relative to Cp T fe 
can "be neglected. In the author's experiments the quanti- 
ties G-/3P and p were measured directly. In the large 
tube preceding the experimental section the inlet tempera- 
ture of the air was measured "by a resistance thermometer. 



♦Tor', the laws of averaging of the separate quantities, see 
the author's work on hydraulic resistance for high velocities. 
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The quantity 0, may he calculated up to any section "by 
the water mass flow rate and the temperature rise of the 
water, 

From the plot of Q and the initial air temperature, 
the quantity so— called "stagnation temperature" T m and 
its change alcng the test section have "been computed 

T m " A ~ + T = T 0 - A (5) 

Cp2g ° Cp 



A solution of equations (2) and (4) gives a quadratic 
equation for the flow temperature 



s 2g8p /? 1\ 2 2gCp /r 1\ 2 

A \G E/ A \G E/ 

To calculate T from equation (6) a graphical solution 
"based on methods is described in the work on flow resist- 
ance. The magnitude of flow velocity is computed from 
the difference between stagnation and actual temperatures. 



w = /^(T m ~T) (7) 



From the experiment friction factors have been cal- 
culated under conditions with heat flow. The Bernoulli 
equation for flow of compressible fluids 



wdw c w 2 dl , . 

= - T dp - 1 (8) 

g 2g D 

together with equation (2) determines the friction factor 
\ for compressible flow. 



i = _ 2D d( P + P" S > (9) 



dLpw 



z 



NACA Technical Memorandum No. 1054 



9 



where 

p 1/Vg = air density 
D tube diameter 

With very low velocities and p w s « p or w « gpV = 

(a acoustic velocity, k adiabatic exponent) — that is, 
with velocities far "below acoustic, equation (9) "becomes 
the ordinary relation for | for flow of incompressible 
fluids. 

The quantity (p + p w 2 ) varied all along the tube 
in the experiment much less than did the quantities p 
and impact pressure p w 2 separately. Therefore in cal- 
culating | from the experiments the curves of (p + p w 2 ) 
against tube length were differentiated. The values of I 
have been calculated for the sections L = 700, 900, 1100, 
1250, and 1360 millimeters. 

In figure 5 the experimental values of \ as a 
function of Peynolds number, calculated on the mean 
stagnation temperature in a given section are shown as 
small circles. The curve fits the Nikuradse equation 



| = 0. 0032 + -2i£21_ ( 

Ee 0.E37 



So it is seen that the actual experiments with heat trans- 
fer occurring confirm the previous conclusion made for 
adiabatic streams — that is, that the Nikuradse relation 
for t is applicable to streams of high velocity* 

In the same sections where the friction factors were 
determined, local heat transfer coefficients from the equa 
tion for heat flow for a tube surface element of diameter 
D and length dL were also computed 

&q, = a (T m - T cm )TTD dL (11) 
where T cm is the absolute temperature of the tube wall. 
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The values of dQ were calculated from experimental 
data by use of equation (l). 

The utility of reckoning heat transfer coefficients, 
not on actual fluid temperature hut on stagnation temper- 
ature, has already been well founded by Vittgeo (refer- 
ence 2), M. F, Shirokov (reference 4), and Jung (refer- 
ence 5). Calculation' of the mean coefficient of heat 
transfer a was .made for tube sections L = 0.7 (L/D=50) 
to L = 1.350M. 

The curves of water and wall temperatures against 
tube length approximated straight lines in the section 
up t o L = 1.250M; therefore for stagnation temperature 
and for Q, , there also is, in this section, approximate 
linear dependency on length according to equations (5) 
and (l). To determine the mean a according to equa- 
tion (11) 



a = i r a dL = x r , *3 ( 12) 

' t J tfDL J T m - T cm 

Substituting 'dQ = q, dL, Tjjj-— Tcm — - T o'*' where ' cc, 

b, i" 0 are constant, and after integration 



a - -3 1 ln Is 



or 



a ^ (I3) 

TTr> L T 



where t = T m - T cm 
given section and 



'is the temperature potential in a 
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is the logarithmic mean temperature difference in the 
tube between two transverse sections* ' [ 

• By use of equation (13) values of ' a have been 
calculated for' all these experiments as' well as the 
dimension-Less parameters Nu and Pe characterizing 
the heat transfer, ' ■ 

^ = ■ • • ( i4) 

Pe = (15) 



Here A is the conductivity of the air. 

The Reynolds— Prandt 1 hydrodynamic theory of heat 
transfer predicates a transfer of heat and momentum from 
the core of the stream to the walls by' a process of tur- 
bulent eddy motion of the same stream particles. If, in 
this theory, the resistance to heat transfer Of the thin 
laminar films for gases in tubes is neglected, the follow- 
ing relation between a and \ . is obtained: 

i 

a - wYOp — ' 



or in dimens ionless parameters, 

5H = I 
Pe 8 



(16) 



As Shirokov has demonstrated,, the hydrodynamic theory 
of heat transfer may be extended to streams . of "high veloc- 
ities if for the transfer , not , only of the heat but also 
of equivalent kinetic energy from the core of the stream to 
the walls is allowed. Then equation (16) is valid if a, 
is calculated on the temperature difference based on stag- 
nation temperature. 

The data of these experiments confirm, this relation 
and check with the "data of Tittgeo (reference 2). 
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In ftgure 6 where the abscissa is Nu/Pe and the 
ordinate is' I the experimental data for mean values 
along the test section are marked with small circles and 
the data of Vi-ttgeo with triangles, the straight line is 
drawn as required by equation (16). It may be seen that 
over almost the whole range, of this investigation, the 
data correlates well with equation (16) except in the 
small region of high values of Pe. 

In figure 7 local values of Nu/Pe and I are 
marked. The solid straight line corresponds to equation 
(16) and the dotted line was fitted to the experimental 
data by the method of least squares. The conclusion re- 
garding the dependency of Nu/Pe on' £ for mean values 
is also confirmed by figure 7. 

Usually for gases, heat transfer experiments are 
reported as Ku against Pe, because £ depends on 
Re . and Ee = Pe/Pr for gases is directly proportional 
to Pe -(the Prandtl criterion Pr depends only on the 
number of atoms in the gas). 

In figur'e 8 Pe values are on the abscissa and Hu 
on- the ordinate and the data based on mean values are 
shown as small circle.s. The values of the physical con- 
stants in !Tu and Pe are calculated on mean gas stag- 
nation temperature. It is to be noted that the data, 
calculated on stagnation temperature and not on limiting 
temperature's as is usually done in this field, fit well 
with the ten Bosch equation for air. 

Nu = 0. 0364 Pe 0 ' 75 ( 17) 

This is shown graphically on the same plot by the solid 
curve. This equation is obtained by substitution in equa- 
tion (16) for £ by the relation of Blasius which fits the 
Nikuradse equation for Re .J* 100,000 and it is confirmed 
by multiple experiments at low velocities if the physical 
properties are based on the wall temperature, 

'This is also shown in figure 9 where the small circles 
show the relation between the local values of Nu and Pe 
calculated on stagnation temperature and the solid curve 
corresponds to the ten Bosch formula. 

It should be stated that these experiments have not 
shown any noticeable influence of temperature on a, and 
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I as has "been reported by Jung on stack gases where 
the temperature variation of the gases was much greater 
than in our experiments. 



CONCLUSIONS 



1. The relations for friction and heat transfer 
for cooling hot air streams flowing at velocities up' 

to acoustic velocity 'through straight water— cooled tubes, 
have been studied experimentally, 

2. The range of variables covered values of Nusselt 
numbers, ■ — Nu from 200 to 500 J Pecle't numbers ---Pe from 

50, 000 to 300,000, and Reynolds numbers— Re from 70, 000 
to 420,000; the air temperature varied from 400° C at the 
inlet to 140° C at the exit, and the velocities reached 
the acoustic value (up to 420M/sec) at the exit end of 
the tube. 

3. It was found that within the limits of our ex- 
periments, the calculation of heat transfer coefficients 
based on stagnation temperatures confirmed the relation, 

— — = — , which is derived from the Reynolds— Prandt 1 
Pe 8 

analogy. This is as well confirmed for mean values of 
Nu and Pe in the tube as for point values of Nu and 
Pe and is "by no means dependent on the (Bairstow) Mach 
number (the ratio of stream velocity to acoustic velocity), 

4. The values of £ for flow with heat transfer fit 
the Nikuradse formula as well as % for streams without 
heat transfer. 



{ = 0. 003 2 + 



0. 231 
Re °' 237 



Translation by N. P. Vakar 
and G. C. Williams. 
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Table I-Continued 
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Figure 2.- Scheme of experimental tube. 
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Figure 4.- Temperature distribution along experimen- 
tal tube. 
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